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angular momentum I, resulting in K values increasing with I within each band, even 
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observed recently for the odd-even staggering in octupole bands of light actinides. 
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1 Introduction 



The properties of nuclear systems with octupole deformations [|T| are of current interest 
due to increasing evidence for the presence of octupole instabilities in various regions of 
the nuclear table [0, ^. Furthermore, some "beat" patterns have been observed recently 
for the odd-even staggering (the relative displacement of the odd levels with respect to the 
positions at which they should have been located according to a fit of the even levels by 
the formula E{I) = AI{I + 1), where / denotes the angular momentum) in octupole bands 
of light actinides based on recent experimental data @, ^, calling for a study of the 
interactions which could give rise to such shapes. 

Various parametrizations of the octupole degrees of freedom p|, ^ 0, |Tl]] already exist, 
being a useful tool for understanding the role of the refiection asymmetry correlations and 
for analyzing the collective properties of such systems. Some important questions in this 
direction are: which are the collective nuclear interactions that correspond to the different 
octupole shapes and how do they determine the structure of the respective energy spectra? 
Physically meaningful answers should be obtained by taking into account the simultaneous 
presence of other collective degrees of freedom, such as the quadrupole ones |]T^|. 



In the present work we address the above problems by examining the interactions that 
generate collective rotations in a system with a simultaneous presence of octupole and 
quadrupole deformations. The basic assumption of our consideration is that the rotational 
motion of such a system can be interpreted in first approximation as the motion of a body 
with a stable quadrupole-octupole shape. In this respect our purpose is to examine how 
the nuclear system behaves under collective rotations if the presence of stable quadrupole- 
octupole deformations is assumed. 

Based on the octahedron point symmetry parametrization of the octupole shape [§, ^, 



rO| , we propose a collective Hamiltonian which incorporates the interactions responsible for 
the rotations associated with the different octupole deformations. In addition we take into 
account the quadrupole degrees of freedom and the appropriate higher order quadrupole- 
octupole interaction. Below it will be shown that such a general model Hamiltonian could 
incorporate the basic properties of a nuclear system rotating under the above assumption. 

Although this assumption seems to be rather strong (since the presence of stable oc- 
tupole deformations in nuclei is not yet a well elucidated problem) we suppose that it could 
give a natural possibility to estimate the extent to which some of the observed nuclear 
octupole bands carry the characteristics of the stable octupole shapes. Generally, the pro- 
posed consideration will provide a direct physical insight into the nuclear collective motion 
as far as the shape of the system and the respective moments of inertia are slightly changed 
under the collective motion. Similarly to the case of the pure quadrupole deformations 
this requirement will be naturally satisfied for the low angular momentum region of the 
spectrum which is, from another perspective, accessible for detailed microscopic analysis, 
the length of this region depending on the particular system. 

Furthermore, we expect that in the higher angular momentum regions the approach 
suggested will outline some general properties of the system and thus will provide a relevant 
guide for respective more detailed studies from both microscopic and phenomenological 
points of view. In particular it will be shown that the model formalism developed in 
the present work proposes a schematic explanation of the recently observed pi "beat" 
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patterns for the odd-even staggering in octupole bands of light actinides based on recent 
experimental data 0. 

In Section 2 of the present work the octupole terms of the Hamiltonian, classified by 
the irreducible representations (irreps) of the octahedron group, will be described, while 
the quadrupole terms and the octupole-quadrupole interaction will be examined in Section 
3, along with the bandhead term of the Hamiltonian and the minimization procedure, 
which is a basic ingredient of the present work. In Section 4 the diagonal parts of the 
Hamiltonian will be analyzed and used for the production of schematic odd-even staggering 
patterns, while the same procedure will be repeated including the non-diagonal parts of the 
Hamiltonian in Section 5. In Section 6 an analysis of the model formalism under some 
restrictions on the permitted values of the angular momentum projection on the body-fixed 
2;-axis will be presented, while further tests of the formalism will be performed in Section 
7. Finally, Section 8 will contain discussion of the present results, while in Section 9 a 
summary of the present results and plans for future work will be given. 

2 Parametrization of the octupole deformation 

Our model formalism is based on the principle that the collective properties of a physical 
system in which octupole correlations take place can be expressed by the following most 
general octupole field in the intrinsic (body-fixed) frame 
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^3 = E 




fj,=—3 




This field can be written in the form 
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(2) 



i=l 



where the quantities 
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F2(3) = -iiv^(y3S + i'3-3) + 4^(^31 + V3_i) 

(with = + + z^) belong to irreducible representations (irreps) of the octahedron 
group O. In particular, the first quantity (Eq. (3) ) belongs to the one- dimensional irrep 
A^^ while the next three quantities (Eqs (4)- (6) ) belong to the three-dimensional irrep Fx 
and the last three quantities (Eqs (7)-(9) ) belong to the three-dimensional irrep ¥2- The 
seven real parameters Cq and e^ij) [r = 1,2; i = 1,2,3), appearing in Eq. (2), determine 
the amplitudes of the various components of the octupole deformation. Their relation to 
the coefficients a^'^ has been given in Ref. 0. 

As we have already seen, the quantities appearing in Eqs (3)-(9), when expressed in 
terms of the cartesian coordinates x, y, and z, contain linear combinations of terms cubic 
in the cartesian variables. These specific linear combinations correspond to various octupole 
shapes (as seen from their expressions in terms of the spherical harmonics), and in addition 
correspond to the above mentioned irreps of the octahedron group. Our proposition is the 
following: 

a) We construct a Hamiltonian using the same cubic terms appearing in Eqs (3)- (9), 
but replacing the cartesian coordinates x, y, z by the angular momentum operators I^, ly, 
Iz (with P = II + ly + ). For example, the term z^ is replaced by . 

b) Due to the fact that the operators Ix-, ly, Iz do not commute, while the coordinates 
X, y, z do commute, when making step a) we symmetrize each cubic expression when we 
write it in terms of the angular momentum operators. For example, the term 2;(x^ — y"^) is 
replaced by /.(Jj - J^) + (/J _ /2)/^. 

c) During the procedure described above, the factors appearing in the denominators 
of Eqs (3)- (9) are replaced by P factors. In the final result we normalize with respect to 
J3, i.e. we multiply the results by P, an operation which is equivalent to the transition to 
a unit sphere, a natural thing to do since we are interested in surface shapes. Indeed, this 
operation is equivalent to the multiphcation of the quantities appearing in Eqs (3)-(9) by 
r^, an action which eliminates their radial dependence, leading to a transition to the unit 
sphere. The transition to the unit sphere is the reason that the standard quadrupole and 
octupole operators are defined to be proportional to r^y2/i and r^K^^ respectively 

d) As a result, we obtain a Hamiltonian (as a function of the angular momentum oper- 
ators Ix, ly, Iz) the terms of which correspond to the same octupole shapes which appear 
in the Hamiltonian of Eqs (2)- (9). 
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e) The terms of the resulting Hamiltonian in addition belong to the same irreps of 
the octahedron group as the terms of the original Hamiltonian, appearing in Eqs (3)- 
(9). In other words, through this procedure we determine the octahedron point symmetry 
properties of the system in angular momentum space. 

Our proposition is similar to the procedure used in Ref. [|l^ for the hexadecapole field. 

The following Hamiltonian is then obtained: 

2 3 

Hoct = HA2 + Y.Y.HF.ii) , (10) 



r=l i=l 



with 



= a2-MiJy + iyh)iz + iz{Liy + iyix)] , (n^ 



Hf.W = ^/ii/.(5/f - 3/2) , (12) 
HF,i2) = \h2{^il - 3/J') , (13) 

^F,(3) = ^/l3(5§ - 3/,?) , (14) 



Hf.w = f2i\[Uil - ID + {il - , (15) 
HF,i2) = f22{LP - il - iJl - 1%) . (16) 

^F.(3) = f23{iyil + iliy + il-iyP) (17) 

The Hamiltonian parameters 02 and f^i (r = 1 ,2; i = 1, 2, 3), appearing in Eqs (11)-(17), 
are formally related to the parameters in Eq. @) as follows 





a2 = eoy^, fii = f2i = ^2{i)\l i^, z = 1,2,3. (1^ 

The non-vanishing matrix elements of the operators of Eqs (|TT])-(|I7p in the states with 
collective angular momentum / are given in the Appendix. 

In the Appendix we remark that the operator Hfi{i) (Eq. ([T^ ) ), which corresponds 
to the term Y^q (see Eq. (4) ), characterized by axial deformation, is the only octupole 
operator possessing diagonal matrix elements. Below it will be seen that this operator is 
of major importance in determining the fine structure of collective bands with octupole 
correlations. This fact does not come as a surprise, since it is well known that the Y30 
(axial) deformation is the leading mode in systems with reflection asymmetric shapes (see 
Ref. [Hi for a relevant review). 



3 Inclusion of quadrupole degrees of freedom 

It is known, however, that the use of the pure octupole field of Eq. (|1|) is not sufficient for 
the description of the collective properties of nuclei exhibiting octupole deformation, since 
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the quadrupole deformation is also present. Therefore one has to consider the octupole 
degrees of freedom together with the quadrupole deformations, and in addition one has 
to deal with their coupling. A general treatment of a combined quadrupole-octupole field 
has been given earlier in the framework of a general collective model for coupled multipole 
surface modes |T^]. 

According to the above considerations it is reasonable to suggest that the most general 
collective Hamiltonian of a system exhibiting octupole deformations should also contain the 
standard (axial) quadrupole rotation part 

Hrot = AP + Ail , (19) 

where A and A' are the inertial parameters. In addition, it is reasonable to try to describe 
the coupling between the quadrupole and the octupole degrees of freedom by the follow- 
ing higher order diagonal quadrupole-octupole interaction term, which corresponds to the 
product F20 ^3 0, 

H,oc = /,oc^(15/| - 14/3/2 ^ 3|J4) ^20) 

since it is known that the axial deformation corresponding to Iso is the leading mode in 
systems with reflection asymmetric shapes 0] and, in an analogous way, the axial deforma- 
tion corresponding to 120 is the leading mode in systems with quadrupole deformations. 
The operator of Eq. (20) is also normalized with respect to the multiplication factor /'^, in 
the same way as the operators of Eqs (11)-(17) are, as described in the previous Section. 
In other words, we use the product P ¥20 I30, in order to ensure that all non-quadrupole 
Hamiltonian terms will be of the same order. 

As a result, the total Hamiltonian of the system can be written in the form 

H = Hbh + Hrot + Hoct + Hqoc , (21) 

where 

Hbh = Hq + fkl-^ , (22) 

is a pure phenomenological part introduced in order to reproduce the bandhead energy in 
the form 

Ebh = Eo + fkK , (23) 

were Eq and ft are free parameters. The bandhead energy Ebh is considered proportional 
to K in analogy to what happens in the five- dimensional quadrupole oscillator model (see 
Appendix 6B of Ref. [0), as well as in the standard rotation-vibration model (see section 
6.5 of Ref. [|14|)- The /C-dependence of Ebh plays an important role in our approach, 
since it provides the correct value of the bandhead angular momentum projection K in the 
variation procedure described below. 

We remark that the Hamiltonian of Eq. (^TJ) is not a rotational invariant in general. 
It does not commute with the total angular momentum operators and, as a result, any 
state with given angular momentum / is energy split with respect to the quantum number 
K. Therefore, the physical relevance of this Hamiltonian depends on the possibility to 
determine in a unique way the angular momentum projection. Our basic assumption is 
that K is not frozen within the sequence of states of the collective rotational band. We 
suggest that for any given angular momentum / the projection K should be determined 
so as to minimize the respective collective energy. The resulting octupole band is then 
the yrast sequence of the energy levels produced by our model Hamiltonian. It should be 



mentioned that a similar procedure has been used in Refs p!5| , p!6| in relation to the A/ = 2 
staggering effect in super deformed nuclei. 
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4 The diagonal parts of the Hamiltonian 



As a first step in testing our Hamiltonian we consider its diagonal part 

n = Hbh + Hrot + -f^orf + Hqoc , (24) 

where the operator H^^^ = Hfi(i) represents the diagonal part of the pure octupole Hamil- 
tonian Hoct of Eq. ( |T0| ) . 

The following diagonal matrix element is then obtained 

Ek{I) = Eo + hK + AI{I+l) + AK^ + h,(^-K^-^-KI{I + l)^ 

+ f,oc^{lhK'' -UK''I{I + 1) + ?>KI\I + If) . (25) 

Following the above mentioned assumption for the angular momentum projection 
we determine the yrast sequence E{I) by minimizing for each value of / the expression 
of Eq. ( p5D as a function of integer K in the range — / < K < I. The obtained energy 
spectrum depends on six model parameters: Eq is a constant contributing to the bandhead 
energy; fk is determined in order to correspond to the correct bandhead energy of the 
octupole band for J = 1, in which case one has K = 1 as well (This is appropriately 
modified in the case of bands starting with higher K. For example in K = 5 bands the 
lowest angular momentum is / = 5, and is fixed in order to provide the correct bandhead 
energy for I = K = A and A are the quadrupole inertial parameters which should in 
general correspond to the known quadrupole shape (axes ratio) of the relevant nucleus; /n 
and fqoc are the parameters of the diagonal octupole (Eq. (|T2|) ) and quadrupole-octupole 
(Eq. (^) ) interactions respectively. In what follows we give fixed values to the first four 
parameters (i?o, /fc, A, A) and vary the last two parameters (/n, fqoc) in order to examine 
the influence of the last two terms of Eq. (25) on the odd-even staggering pattern. 

A schematic energy spectrum of this kind, obtained for a sample set of parameter values, 
is given in Table 1. It is seen that the "yrast" values of the quantum number resulting 
from the minimization procedure described above, gradually increase with the increase of 
the angular momentum /. We remark that these values of K correspond to the local minima 
of the energy expression of Eq. (|25|) as a function of K, which is illustrated in Fig. 1 for 
the same set of parameters. We see that these minima appear for positive values of K, are 
well determined, and their depth increases with the increase of the angular momentum /. 

These results indicate that the rotational motion of a system with stable quadropole- 
octupole deformation is associated with complex angular momentum dynamics, due to 
the complex shape contributions to collectivity. More precisely, the following theoretical 
proposition is formed: For this kind of rotational motion the increase in the total angular 
momentum of the system is associated with an energetically favorable increase in its third 
projection K on the body fixed axis. This means that the vector of the angular momentum 
deviates "step by step" from the "x-y" body-fixed plane (which is perpendicular to the body- 
fixed "z"-axis) as its magnitude increases. In other words, the higher angular momenta will 
be attended by stronger precession, or "wobbling motion" of the system. These findings 
resemble the recent results [ITTI obtained within the framework of a self-consistent cranked 
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HFB approach. Application of the method to some Os isotopes has demonstrated that 



high-X bands become important in the high angular momentum region [[I^ . 

Such a behavior of the spectrum could also be interpreted as a multiband-crossing 
phenomenon, since the obtained yrast sequence can be considered as the envelope of the 
curves with different values of the quantum number K, as it is illustrated in Fig. 2. This 
interpretation can, however, be considered as a mathematical description, not necessarily 
implying the physical existence of several bands with different values of K. 

Thus, our schematic consideration suggests that the rotational motion of a system with 
stable quadrupole-octupole deformations invokes i^-values essentially higher than the ones 
usually considered in microscopic studies. We remark that the situation considered here is 
essentially different from the case of pure quadrupole rotations, where the bandhead value 
always provides the minimum of the collective rotational energy along the whole band (for 
example Erot = AI[I + 1) + A' K"^) . Thus the obtained wobbling motion appears as an effect 
of the higher multipolarity of the considered collective interactions. The extent to which 
the K-values appearing in our procedure provide physically reasonable interpretations of 
various nuclear collective modes will be considered in Section 6, where the limitation to 
K <3 will be studied. 

Now we can examine the fine structure of the collective bands obtained in the above 
model procedure. We see in Table 1 that the iC-values characterizing the levels with 
increasing / are grouped in couples containing two successive values of I (one even and 
one odd). This fact implies the presence of an odd-even staggering effect, as we shall 
immediately see. A measure of the odd-even staggering effect is the quantity 

Stg{I) = 6AE{I) - AAE{I - 1) - AAE{I + 1) + AE{I + 2) + AE{I - 2) , (26) 

with 

AE{I) = E{I + 1) - E{I). (27) 

The quantity Stg{I) is proportional to the discrete approximation of the fourth derivative 
of the function AE{I), i.e. proportional to the fifth derivative of the energy E{I), and 
is able to demonstrate fine deviations from the pure rotational behavior, as it has been 
demonstrated in Ref. for the octupole bands of several light actinides. An analog of 
this quantity has been introduced earlier |jl8|, |19|, ^ for the study of AI = 2 staggering in 



nuclear superdeformed bands. 

The way in which the appearance of the K-values in couples for successive values of / 
leads to odd-even staggering can be seen from Table 1, using the second term of Eq. (25) 
as an example. In Table 1 we see that couples of K values appear for / > 8. In this region 
the contribution of the second term of Eq. (25) to AE{I) is zero for even values of /, while 
it is for odd values of /. Then we easily see from Eq. (26) that the contribution of the 
second term of Eq. (25) to Stg{I) is +8/^ for odd values of /, while it is — 8/^ for even 
values of /. Therefore the second term of Eq. (25) leads to odd-even staggering of constant 
amplitude. The 4th, 5th and 6th terms of Eq. (25) depend nonlinearly on K and/or J, and 
therefore give odd-even staggering of varying amplitude. In the case of the 4th term, it is 
easy to see that it gives staggering with amplitude increasing linearly as a function of I. 
(The terms Eq and AI{I + 1) do not contribute to the odd-even staggering, as it is known 
from Ref. [§.) 
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Several examples of odd-even staggering occurring from Eq. (26) for different sets of 
parameter values used in Eq. (25) are given in Fig. 3. The parameters Eq and fk have 
been kept constant in all parts of Fig. 3, while the parameters A and A' have been kept 
constant in all parts of Fig. 3 except the last one (Fig. 3(f) ). In contrast, the parameters 
/ii and fgoc have been given different values in the various parts of Fig. 3. 

Fig. 3(a) illustrates a long odd-even staggering pattern which looks similar to the 
"beats" observed in the octupole bands of some hght actinides § (with ^^^Ra, ^^^Ra and 
^^^Ra being probably the best examples), a difference being that in the realistic cases of 
Ref. 1^ the amplitude of the "beats" seems to be decreasing with increasing /, while in 
Fig. 3(a) the opposite seems to hold. 

In Fig. 3(b) the increased values of /n and fgoc provide a wide angular momentum 
region (up to / ~ 40) with a regular staggering pattern. 

The further increase of fgoc results in a staggering pattern with different amplitudes, 
shown in Fig. 3(c). 

The further increase of fu and fgoc leads to a staggering pattern with many "beats", 
as shown in Fig 3(d). Notice that in Fig. 3(d) the first three "beats" are completed by 
/ ^ 40, while in Fig. 3(a) the first three "beats" are completed by / ~ 70. 

Fig. 3(e) illustrates what happens in the case of vanishing quadrupole-octupole interac- 
tion term {fgoc = 0), keeping the rest of the parameters the same as in Figs 3(b) and 3(c). 
We see that in the present case the "beat" effect occurs very frequently, while in Figs 3(b) 
and 3(c) the change was much slower. 

Finally, an example with almost constant staggering amplitude is shown in Fig. 3(f). 
It resembles the form of the odd-even staggering predicted in the SU(3) limit of various 
algebraic models (see Ref. for details and relevant references). It also resembles the 
odd-even staggering seen in some octupole bands of light actinides (with ^^^Rn being 
probably the best example). 

5 The non-diagonal parts of the Hamiltonian 

After seeing in the previous section the main features of the diagonal part of the Hamil- 
tonian, we can now focus our attention on the general Hamiltonian of Eq. ([21|), including 
the various non-diagonal terms given in Eqs (|n|), (|TE|)-(|T^. The main problem in this 
case is the fact that K is in general not a good quantum number. We can, however, make 
analysis for small values of the respective parameters (02, /12, /13, /21, /22, /23), which 
keep K "asymptotically" good. Keeping these parameters small means that we use a weak 
ii'-bandmixing interaction, which guarantees that for any explicit energy minimum appear- 
ing in the diagonal case the corresponding eigenvalue of the perturbed Hamiltonian will be 
uniquely determined. Thus we are able to obtain for the perturbed Hamiltonian a K-mixed 
yrast energy sequence analogous to the one we got in the previous section for the diagonal 
part of the Hamiltonian. 

Our numerical analysis of the Hamiltonian eigenvector systems shows that the param- 
eters of the non-diagonal terms should be smaller by an order of magnitude in comparison 
to the value of the parameter fu. In addition, by examining the corresponding matrix 
elements in the Appendix, we deduce that the following pairs of non-diagonal terms give 
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identical contributions to the energy spectrum: and Hp^i^iy, i^Fi(2) and -^^1(3); -^^^2(2) 

and Hf2{3)- Therefore it suffices to keep from now on only the terms -^^1(2), -^^2(1)) -^F2{2)- 
In Fig. 4 two staggering patterns in the presence of X-bandmixing terms are illustrated. 
In the calculation we have included, as we have already mentioned, the non-diagonal terms 
Hp^(2)-, Hp^(^ij, and Hf2{2)i along with the already considered diagonal Hamiltonian of Eq. 
p^). The three non-diagonal terms have been included with two different sets of parameters 
(both obeying the above mentioned condition that the values of the parameters accompa- 
nying the non-diagonal terms should be at least an order of magnitude smaller than the 
value of the parameter /n), shown in Figs 4(a) and 4(b), while the parameters of the di- 
agonal part have been kept the same as in Fig. 3(b) (and in Table 1). Comparing Figs 
3(b) and 4(a) we see that the non-diagonal terms affect more severely the higher angular 
momentum region by decreasing the staggering amplitude. Moreover, the larger values for 
the parameters of the non-diagonal terms, shown in Fig. 4(b), reduce more seriously the 
staggering pattern in the higher angular momentum region, as one can see by comparing 
Figs 3(b) and 4(b) and noticing the different scales on the vertical axes of the two figures. 
The pattern appearing in Fig. 4(b) resembles the experimental situation in ^^^Rn and ^^^Th 
P (odd-even staggering with amplitude decreasing as a function of J). 

6 \jOW-K {K < 3) analysis of the model formalism 

A discussion on the X-values appearing in our model consideration is appropriate in this 
place. 

Usually from microscopic point of view the eigenvalues of the total angular momentum 
projection 1^ (which, in case of axial symmetry, coincides with the intrinsic momentum 
projection J^) are restricted to K =< 3. This requirement is well justified in the case of 
octupole vibrations of quadrupole deformed nuclei. It reflects the shell model concept that 
the octupole degrees of freedom could be generated by an octupole-octupole A = 3 inter- 
action (where A stands, as usual, for the multipolarity) which couples single-particle states 
of opposite parity [Q. In particular, this assumption gives a satisfactory interpretation for 
the fine structure of negative parity rotational bands built on octupole vibrations. 

In view of the above, it is interesting to examine our formalism in the case when the 
space of the i^-values is restricted to Kmax = 3. Indeed such a test is important especially 
for the low angular momentum regions, where the intrinsic structure of octupole degrees of 
freedom is microscopically well studied. 

Along these lines, we performed a schematic calculation including the values K = 0, 1, 
2, 3 in the angular momentum region up to / = 10, as well as another calculation including 
the values = 0, 1, 2, 3 in the angular momentum region up to J = 16. The yrast 
sequences obtained are given in Tables 2 and 3 respectively, while the resulting staggering 
patterns are shown in Fig. 5. We see that our formalism provides a regular staggering 
pattern for the states up to / = 10 [Fig. 5(a)] and a "beat" pattern for the spectrum up 
to / = 16 [Fig. 5(b)]. These examples (obtained for the parameter values given in Tables 2 
and 3, respectively) illustrate that our approach is capable to provide reasonable staggering 
patterns under the restriction to < 3, at least in the low angular momentum region 
with / = 10-16. In other words, the schematic model gives reasonable results under the 
restriction to K < 3, at least in the microscopically accessible angular momentum regions. 
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It is, however, important to remark that our approach suggests a rather more extended 
treatment of the collective octupole degrees of freedom, which is expected to be useful 
in view of the increasing bulk of data indicating possible stable octupole deformations in 
several nuclei, the light actinides providing the best examples. Some of these bands (the 
octupole bands in ^^"^"^^^Ra and 224-228rp]-^^ example) possess rather well pronounced 
rotational structures (see Table 1 in Ref. [^). Our proposition is that in these cases the 
concept of a rotating quadrupole-octupole shape should be more evident. In these cases, 
in contrast to the pure octupole vibrations, it is not a priori clear why the combined 
quadrupole-octupole intrinsic configurations should be restricted to K < 3 projection val- 
ues. Namely for these kinds of rotational motion our collective formalism suggests that K 
could be higher than 3. One can easily see that for the relatively well deformed nuclei, such 
as ^^^~^^^Ra and ^^le Nilsson orbitals do not exclude the presence of reasonable 

K > 3 intrinsic configurations This remark is of particular importance for the higher 
angular momentum region (/ > 10-16), where pair breaking effects are possible. 

Another important question to be examined is whether the energy gain in case of yrast 
rotations with K > 3 is larger than the intrinsic energy necessary to raise K above 3. There 
is no a priori answer to this question, which would be an interesting problem for future 
combined efforts from both collective and microscopic points of view. Here we just remark 
that for the collective motion considered the energy gain rapidly increases with the increase 
of the angular momentum (see the depth of the minima in the curves of Fig. 1). 

7 Further tests of the model formahsm 

The restriction to K < 3 values, used in the previous section, demonstrated that odd-even 
staggering and "beat" patterns can indeed appear in this special case. This particular 
collection of K values {K = 0, 1, 2, 3) has certain probability to appear low in energy 
in the case of nuclei exhibiting octupole deformation, but other possibilities also exist. In 
fact, it is well known from experiment that the various K values are built up through 
the pair-breaking mechanism. As a result, K is not a smoothly increasing function (while 
/ is). Which K values will appear lower in energy in each nucleus depends strongly on 
the microscopic structure of the nucleus. One is therefore confronted with the following 
questions: 

1) Given a collection of K values, which appear low in energy in a specific nucleus 
because of microscopic reasons, is the present formalism predicting odd-even staggering or 
not? 

2) What is the meaning of minimization with respect to K in this case? 

3) In which order do the various allowed values of K appear as I is increasing? 

We have performed several numerical tests, which need not to be reproduced here in 
detail, in order to answer these questions. The following conclusions have been reached: 

1) Given a collection of K values which are non- successive integers (for example, K = 5, 
7, 9, 10, 15), odd-even staggering and "beat" patterns do appear for certain values of the 
Hamiltonian parameters. One is therefore persuaded that in principle the model formalism 
can produce "beat" patterns for any microscopically directed collection of K values. 

2) In the just mentioned case, K is not a continuous variable. However, one can still 
minimize the energy with respect to K, in the following sense: Given a collection of K 
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values, for each value of / one determines (among the members of the collection) the K 
value for which the energy is minimum. 

3) In all cases the allowed values of K appear in increasing order as / is increasing, 
similarly to what occurs in Tables 1-3. 

In all of the above we have been assuming that K takes only integer values, i.e. that K 
is a good quantum number. This is an assumption of the model, which is quite usual in the 
realm of phenomenological models (but not in the case of microscopic models), in which K 
is assumed to be a good (or asymptotically good) quantum number. One might wonder, 
however, how strongly the appearance of odd-even staggering and "beat" patterns in the 
present model depends on this assumption. In order to confront this question, we have 
performed the following tests: We repeated the calculations concerning Figs. 3(a) and 3(b) 
allowing K to vary not with a step of 1.0, but with a step of 0.5, or 0.1, or 0.01, or 0.0001, 
as well as allowing i^' to be a continuous variable. The results need not be reproduced 
here in full detail, but the conclusion was that in both cases the odd-even staggering drops 
dramatically as the step of the variation of K decreases, and practically vanishes when K is 
a continuous variable. One therefore concludes that the appearence of odd-even staggering 
and "beat" patterns in the present model strongly depends on the assumption that K is 
a good (or asymptotically good) quantum number, taking integer values, an assumption 
which is rather common among phenomenological models. In the other hand, as it has 
already been mentioned in the beginning of this section, practically any microscopically 
directed set of integer values of K can lead to the appearance of odd-even staggering and 
"beat" patterns within the framework of the present model. 

In conclusion, the assumption that X is a good (or asymptotically good) quantum num- 
ber, therefore taking integer values only, is a fundamental one as far as the appearance of 
odd-even staggering and "beat" patterns within the present model is concerned. Once K 
is assuming only integer values, odd-even staggering and "beat" patterns can in princi- 
ple be constructed (for appropriate parameter values in the Hamiltonian) for any set of 
microscopically directed K values lying low in energy. 

8 Discussion 

The staggering patterns illustrated so far (Figs 3 and 4) cover almost all odd-even staggering 
patterns seen in nuclear octupole bands. The amplitudes obtained for the sets of parameters 
considered vary up to 300 keV. Staggering patterns with larger values of Stg{I) can be easily 
obtained for different parameter sets. On this basis it is reasonable to assume that the model 
parameters can be adjusted appropriately so as to reproduce the staggering patterns seen in 
the octupole bands of the light actinides |^ , as well as in rotational negative parity bands 
with K > 1 built on octupole vibrations. In the fitting procedure one should take into 
account the minimization of the energy with respect to K for each given value of /. Work 
in this direction is in progress. 

At this point the following comments on the structure of the collective interactions used 
and the related symmetries are in place: 

1) The above mentioned fact that the six non-diagonal terms of the octupole Hamilto- 
nian can be arranged into three pairs, the terms belonging to the same pair giving equal 
contributions, indicates that only four terms of the octupole Hamiltonian (the diagonal 
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term plus three appropriately chosen non-diagonal terms) suffice for the determination of 
the energy spectrum. This result reflects the fact that in the intrinsic frame of reference 
three octupole degrees of freedom, from the total of seven ones, are related to the orienta- 
tion angles. In the present case we chose in Sec. 5 and in Fig. 4 to keep the diagonal term 
Hfi{i) and the non-diagonal terms Hfi{2), -^F2(i) and Hp,^(2). Our analysis (related to the 
collective rotations of the system) gives a natural way for determining the four collective 
octupole interaction terms which give independent contributions. 

2) The symmetry of the various combinations of spherical harmonics appearing in Eqs 



(3)-(9) has been considered in detail in Ref. [IC]. The I30 term is axially symmetric, i.e. it 
has the symmetry of the Doo group, while the term Y^i — Y^^i has the symmetry of the 
group, the term 132 + ^3-2 has the symmetry of the T^ group, and the term Fss — 13-3 has 
the symmetry of the group. /^From the correspondence between Eqs (3)- (9) and Eqs 
(11)-(17) we see that the diagonal term Hp^ii) is axially symmetric, while the non-diagonal 
terms Hfi{2)-, Hf^(t^) and Hf2{2) (mentioned in the previous item 1) ) are constructed by 
using the combinations (Isi - Ys-i), (^32 + ^3-2), and (Y33 - Is-s)- 

3) Concerning the role of the various terms in the odd-even staggering effect, our 
schematic results of Sees. 4 and 5 indicate that the non-diagonal K-bandmixing inter- 
actions suppress the staggering pattern, while the axially symmetric term Hfi(i) is able to 
provide a "beat" staggering behavior (see Fig. 3(e)). The quadrupole-octupole term Hgoc 
does influence the staggering pattern, providing wider angular momentum regions with 
regular staggering. In short, the axially symmetric term seems to be the most important 
one (among the terms involving the octupole degrees of freedom) for the production of 
"beat" behavior in the odd-even staggering pattern. However, it should be remembered 
at this point that the octupole and octupole-quadrupole terms are not the only ones con- 
tributing to the odd-even staggering effect. As we have already mentioned in Sec. 4, the 
term fkK, coming from the bandhead energy, makes to the odd-even staggering a contribu- 
tion of constant amplitude, while the A'K"^ term, coming from the quadrupole part of the 
Hamiltonian, makes a contribution to the odd-even staggering with an amplitude which is 
increasing linearly as a function of /. It should be noticed, though, that these remarks are 
based only on a few schematic calculations and in no way are final conclusions. In order 
to reach safer conclusions about the relative importance of the various terms in giving rise 
to a "beat" behavior of the odd-even staggering, one has to perform detailed ffis to the 
octupole bands of the light actinides, as mentioned above. 

4) As has already been commented is Sees. 4 and 6 the most general schematic results 
we have obtained suggest that in the high angular momentum region some high K values 
should be involved (as a result of the minimization of the energy with respect to K for each 
given value of /, which is an important ingredient of our approach). The fact that the K 
quantum number is not a good quantum number (not even approximately) of the relevant 



states has been realized long ago |2^. In microscopic calculations in the rare earth 



region (with 152 < A < 190) [^, in which the values K = 0, 1, 2, 3 have been included, 
it has been seen that in the beginning of the region the values K = 0, 1 are important for 
the lowest 3~ state, while in the middle of the region the values = 1, 2 are important 



and in the far end of the region the values i^' = 2, 3 are important [22|. These results are 



consistent with our low-i^ model analysis given in Section 6. The same authors have dealt 
with the actinide region {A > 222) in Ref. p3|. One of the authors of Refs p2| , p3| in Ref. 
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23] finds that the restriction to < 3 is not justifiable for large energies. These findings 



are in agreement with our results given in Table 1 as well as with our comments in the end 
of Sec. 6. 

5) Concerning the analysis in Sec. 5 and the above comments 3) and 4), we remark that 
the restrictions imposed on the non-diagonal Hamiltonian terms (keeping K asymptotically 
good) are rather strong and refiect the particular physical assumptions of our consideration. 
However, for a more general quantum mechanical system there is no principal reason to 
restrict the analysis by assuming the presence of small non- axial deformations only |25|. In 
this direction an extension of the present work not limited to small non-diagonal contribu- 
tions could be of interest. 



9 Conclusions and outlook 

In summary, we have considered a Hamiltonian involving octupole and quadrupole terms, 
along with an octupole-quadrupole interaction and a i^"-dependent bandhead term. The 
octupole terms have been classified by using the irreducible representations (irreps) of the 
octahedron group O. By minimizing the energy with respect to the angular momentum pro- 
jection K for each given value of the angular momentum I we have reached the conclusion 
that K is increasing with increasing /, even in the case in which K is allowed to assume only 
a resticted set of microscopically dictated values, a result which can be interpreted as cor- 
responding to a wobbling motion. Various terms of the Hamiltonian give rise to odd-even 
staggering depending on / in different ways, in this way making the Hamiltonian able to 
produce several different odd-even staggering patterns, some of which have been observed 
in the octupole bands of light actinides. In order to examine the relative importance of the 
various terms giving rise to odd-even staggering, detailed fits of the octupole bands of the 
light actinides should be performed, a procedure which is not very simple because of the 
above mentioned minimization procedure involved in it. In addition, the applicability of the 
present formalism to rotational bands with K > 1 built on octupole vibrations, which also 
demonstrate odd-even staggering effects, should be examined. The present Hamiltonian, 
because of the variety of odd-even staggering patterns it can produce, seems to be a good 
starting point for systematizing the different odd-even staggering patterns seen in octupole 
bands, as well as in rotational bands built on octupole vibrations. Work in this direction 
is in progress. 
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Appendix 

Non-zero matrix elements of tlie operators of Eqs (p!TD -([l7|) in the states \IK) with total 
angular momentum / and projection K, with X = I {I + 1): 



{iK + 2\Ha, 


\IK 


l^IK - 


2\Ha, 


\IK 


(iK 




\IK 


IK+1 


Hf^{2) 


\IK 


IK-1 




\IK 


IK + 1 


HFiiS) 


I\K 


IK-1 




\IK 


IK + 2 


Hf2{i) 


\IK 


IK -2 


Hf2{i) 


\IK 


IK + 1 


Hf2{2) 


\IK 


IK-1 


Hf2{2) 


\IK 


IK+1 


Hf2{3) 


\IK 


IK-1 


Hf2{3) 


\IK 



a2{K +1)JX-{K + 1){K + 2)JX - K{K + 1), 



--a,iK -1)^X-{K- 1){K - 2)^X - K{K - 1), 
i/ni^(5K'-3X), 



-/i2(3X - 15^2 - IbK - 10)^X - K{K + 1), 



— /i2(3X - IbK^ + 15K - 10)^X - K{K - 1), 
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/i3(3X - 15K' - 15K - 10) J X - K{K + 1) 



-— /i3(3X - 15^2 + 15K - 10) JX - K{K - 1), 
io * 



-f2i{K + 1)^X -{K+1){K + 2)^X - K{K + 1), 
- 1)^X -{K-1){K- 2)^X - K{K - 1), 



-MX - 5K^ -5K-2)^X- K{K + 1), 
1 



-/22(X - 5K' + 5K- 2)JX - K{K - 1) 



--/23(X - 5K' -5K-2)JX- K{K + 1), 



-/23(X - 5K' + 5K- 2)JX - K{K - 1). 
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Table 1: The "yrast" energy levels, E{I) (in keV), and the respective values (in h) 
obtained from Eq. ( 125|) for the parameter set Eq = 500 keV, fk = —7.5 keV, A = 12 keV, 
A' = 6.6 keV, /n = 0.56 keV, fgoc = 0.085 keV, by minimizing the energy with respect to 
K for each given value of /. See Sec. 4 for further discussion. 



/ 


E{I) 


K 


I 


E{I) 


K 


/ 


E{I) 


K 


1 


522.772 


1 


13 


2335.81 


5 


25 


5453.12 


11 


2 


568.327 


1 


14 


2576.57 


6 


26 


5694.49 


12 


3 


637.095 


1 


15 


2827.57 


6 


27 


5935.50 


12 


4 


728.710 


1 


16 


3082.36 


7 


28 


6157.50 


13 


5 


840.857 


2 


17 


3344.94 


7 


29 


6378.29 


13 


6 


971.155 


2 


18 


3608.18 


8 


30 


6575.37 


14 


7 


1123.22 


2 


19 


3877.05 


8 


31 


6770.62 


14 


8 


1288.09 


3 


20 


4143.16 


9 


32 


6937.23 


15 


9 


1472.71 


3 


21 


4413.03 


9 


33 


7101.62 


15 


10 


1668.56 


4 


22 


4676.45 


10 


34 


7232.21 


16 


11 


1880.56 


4 


23 


4942.01 


10 


35 


7360.44 


16 


12 


2101.68 


5 


24 


5197.18 


11 


36 


7449.45 


17 
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Table 2: Same as Table 1, but for the parameter set Eq = keV, fk = keV, A = 12 keV, 
A' = 9.96 keV, fu = 0.75 keV, fgoc = 0.15 keV. The restriction K <3has been imposed. 
See Sec. 6 for further discussion. 



/ 


E(I) 


A' 


/ 


E(I) 


K 


1 


24.00 





6 


487.94 


1 


2 


72.00 





7 


640.19 


2 


3 


144.00 





8 


811.21 


2 


4 


237.97 


1 


9 


994.40 


3 


5 


351.64 


1 


10 


1194.18 


3 



Table 3: Same as Table 1, but for the parameter set Eq — 100 keV, fk — 0.896 keV, 
A = 19.2 keV, A' = 18.24 keV, fu = 0.49 keV, fgoc = 0.074 keV. The restriction K <3 
has been imposed. See Sec. 6 for further discussion. 



/ 


m 


K 


/ 


m 


K 


1 


138.40 





9 


1803.42 


1 


2 


215.20 





10 


2177.43 


1 


3 


330.40 





11 


2580.53 


2 


4 


484.00 





12 


3017.25 


2 


5 


676.00 





13 


3490.32 


2 


6 


905.66 


1 


14 


3991.13 


3 


7 


1167.54 


1 


15 


4521.73 


3 


8 


1466.79 


1 


16 


5087.69 


3 
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Figure Captions 



Figure 1. The diagonal energy matrix element Ek{I) (in MeV), Eq. (^), is plotted as 
a function of K (in K) for / = 1, 2, 10, for the parameter set Eq = 500 keV, fk = —7.5 
keV, A = 12 keV, A' = 6.6 keV, fu = 0.56 keV, fgoc = 0.085 keV. 

Figure 2. The diagonal energy matrix element Ek{I) (in MeV), Eq. (^), is plotted as a 
function of / (in %) for K = 10, 11, 12, 13, for the parameter set of Figure 1. 

Figure 3. Odd-even staggering patterns (calculated using Eq. ( p6D ) as functions of /, 
obtained from the diagonal Hamiltonian of Eq. (|^) for several different sets of model 
parameters, given on the figures. Part (b) corresponds to the same set of parameters as 
Table 1 and Figs 1 and 2. The first four parameters remain the same in all parts of the 
figure, except (f). 

Figure 4. Odd-even staggering patterns (calculated using Eq. (^) ) obtained by adding 

the three non-diagonal terms Hfi{2) (Eq- (|T3| ) ), -^F2(i) (Eq. (|T5D ), and Hf2{2) (Eq. (|16D ) 
to the diagonal Hamiltonian of Eq. (p^, for two different sets of model parameters, given 
on the figures. The first six parameters are the same as in Figs 1, 2 and 3(b). 

Figure 5. Same as Fig. 3, but with the restriction K < 3 imposed in the angular 
momentum regions: (a) / < 10 and (b) / < 16. The parameter sets correspond to Tables 
2 and 3 respectively. 
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